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ABSTRACT 

Wc introduce a new descriptor of the weblike pattern in the distribution of galaxies and matter: the 
scale dependent Betti numbers which formalize the topological information content of the cosmic mass 
distribution. While the Betti numbers do not fully quantify topology, they extend the information 
beyond conventional cosmological studies of topology in terms of genus and Eulcr characteristic used 
in earlier analyses of cosmological models. The richer information content of Betti numbers goes 
along with the availability of fast algorithms to compute them. When measured as a function of scale 
they provide a "Betti signature" for a point distribution that is a sensitive yet robust discriminator 
of structure. The signature is highly effective in revealing differences in structure arising in different 
cosmological models, and is exploited towards distinguishing between different dark energy models 
and may likewise be used to trace primordial non-Gaussianities. 

In this study we demonstrate the potential of Betti numbers by studying their behaviour in simu- 
lations of cosmologies differing in the nature of their dark energy. 

Subject headings: Cosmology: theory — large-scale structure of Universe — Methods: data analysis 
— Methods: numerical 



1. INTRODUCTION 

The large scale distribution of matter mapped by 
galaxy surveys like the 2dF , SPSS and 2MASS red- 
shift surveys (iColless et al.l [200l iGott et al.l [20051 : 
iHuchra et aLll2005| ) shows a complex network of inter- 
connected filamentary galaxy associations. This nct- 
work, which has become known as the Cosmic Web 
(|Bond et al.lll996l ). contains structures from a few mega- 
parsecs up to tens and even hundreds of megaparsccs in 
size. Galaxies and mass exist in a wispy web-like spa- 
tial arrangement consisting of dense compact clusters, 
elongated filaments, and sheet-like walls, amidst large 
near-empty voids, with similar patterns existing at ear- 
lier epochs, albeit over smaller scales (for a review, see 
Ivan de Wevgaert fc BoiIdl[20M ). 

Given the structural richness of the Cosmic Web, in 
this paper we address the question of how we can exploit 
the information content in the structure of the Cosmic 
Web to differentiate cosmological models having differ- 
ent dark energy content. Different dark energy models. 
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viewed at the same redshift, show subtly different struc- 
tures in the matter distribution. This is simply a conse- 
quence of the different growth rates for cosmic structure 
in the different models. 

Despite the multitude of descriptions, it has remained 
a major challenge to characterize the structure, geometry 
and topology of the Cosmic Web. Many of the attempts 
to describe, let alone identify, the features and compo- 
nents of the Cosmic Web have been of a rather heuristic 
nature. Moreover, the overwhelming complexity of both 
the individual structures as well as their connectivity, the 
lack of structural symmetries, its intrinsic multi-scale na- 
ture and the wide range of densities in the cosmic matter 
distribution has prevented the use of simple and straight- 
forward methods. 

Here we introduce a new topological measure that 
is particularly suited to differentiating web-like struc- 
tures. To this end, we advance the topological charac- 
terization of the Cosmic Web to a more complete de- 
scription, the homology of the distribution as measured 
by the scale-dependent Betti n umbers of the sample 
(liildelsbrunner fc Mueckd 119941: lEdelsbrunner fc Harerl 
120101: IRobinsI 120061 Ivan de Wevgaert et al.ll20ia I2011D . 
The Betti numbers underlie genus analysis and may 
be regarded as fundamental topological structure in- 
dicators that specifically characterise the homology of 
the distribution. While a full quantitative character- 
ization of the topology of the cosmic mass distribu- 
tion may not be feasible, the homology is an attrac- 
tive compromise, providing a usefully detailed summary 
measurement of topology with relative ly low compu- 
tational overhead (e.g. [Pclfinado fc Edelsbrunnerl 119931 : 
lEdelsbrunner fc Mueckdll994f) . 

1.1. Dark energy and the Cosmic Web 

The parameters for what might now be called "the 
standard model for cosmology" have been established 



2 




Fig. 1. — Slices of thickness 0.2h ^Mpc from two cosmological 
models: ACDM (left), SUGRA (right). Top: 200/1-^ x200h-iMpc 
slices. Middle & Bottom: zoom-in, 40h~^ x40/i"^Mpc slices. Mid- 
dle: density field. Bottom: halo distribution, superimposed. 



with remarkable precision. However, there remains the 
great mystery of the nature of the so-called "dark en- 
ergy" (DE) which appears to make up some 73% of the 
total cosmic energy density. The simplest model for the 
dark energy is Einstein's cosmological constant, A: this 
makes a constant, time independent, contribution to the 
total energy density in the Friedman-Lemaitre equations. 
Models based on this are referred to as "ACDM models" . 
However, there are numerous, possibly more plausible, 
models in which the dark energy evolves as a function of 
time. These models are generally described in terms of a 
time or redshift dependent function w{z) that describes 
the history of equation of state of the dark energy com- 
ponent: 

p = w{z)pc^ (1) 

Different models for dark energy produce different func- 
tions w{z). The Einstein cosmological constant corre- 
sponds to w{z) = —1. 

The aim of this paper is to compare, at redshift z = 
0, the topology of the cosmic structure in simulations 
of three cosmological models having different dark en- 
ergy content that evolved from the same initial condi- 
tions (fig. [ij. The equation of state of the dark en- 
ergy models are those of the standard ACDM model 



and two quintessence models. The latter assume that 
the Universe contains an evolving quintessence scalar 
field, whose energy content manifests itself as dark en- 
ergy. The two quintessence models are the Ratra-Peebles 
(RP) model a nd a SUGRA model ( f or a de tailed de- 
scription. see iRatra fc PeeblesI 119881 lAinendola .20001: 



Brax fc MartinlllQQa iPerrotta et al.ll2000l: Ide Boni et alj 
2OIIHB0S et al.ll201lD . 



2. TOPOLOGY AND HOMOLOGY 
2.1. Beta Numbers 

Formally, homology groups and Betti numbers char- 
acterize the topology of a space in terms of the 
relationship betwee n the cycles and boundaries we 
find in the space^" fMunkrej Il995t iZomorodianI 120051 : 
lEdelsbrunner fc Har er 2010). 

For example, if the space is a d-dimensional manifold, 
M, we have cycles and boundaries of dimension p from 
to d. Correspondingly, M has one homology group 
iJp(M) for each of d -I- 1 dimensions, Q < p < d. By 
taking into account that two cycles should be considered 
identical if they differ by a boundary, one ends up with a 
group iJp(M) whose elements are the equivalence classes 
of p-cycles. The rank of the homology group Hp{M) is 
the p-ih Betti number, /3p = /3p(M). 

In heuristic - and practical - terms, the Betti numbers 
count topological features and can be considered as the 
number of p-dimensional holes. When talking about a 
surface in 3-dimensional space, the zeroth Betti number, 
/3o, counts the components or "islands", the first Betti 
number, counts the tunnels, and the second Betti 
number, /32, counts the enclosed voids. All other Betti 
numbers are zero. 



2.2. Genus and the Euler characteristic 

Numerous cosmological studies have considered 
the genus of the isodensity surf aces defined by 
the megaparsec galaxy distribution ( Gott et all 119861 : 
iHamilton eFall Il986t iChoi et all 120101) . which specifies 
the number of handles defining the surface^^. The genus 
has a simple relation to the Euler characteristic, x, of 
the isodensity surface. Consider a 3-manifold subset M 
of the Universe and its boundary, 9M, which is a 2- 
manifold. With 9M consisting of c = /3o(9M) compo- 
nents, the Gauss-Bonnet Theorem states that the genus 
of the surface is given by 



G-c-ix(5M), 



(2) 



Assuming the space is given as a simplicial complex, a p-cycle 
is a p-chain with empty boundary, where a p-chain, 7, is a sum of 
p-simplices. The standard notation is 7 = ^aiai, where the ai 
are the p-simplices and the ai are the coefficients. For example, 
a 1-cyclc is a closed loop of edges, or a finite union of such loops, 
and a 2-cycle is a closed surface, or a finite union of such surfaces. 
Adding two p-cycles, we get another p-cycle, and similar for the 
p-boundaries. Hence, we have a group of p-cyclcs and a group of 
p-boundaries. 

For consistency, it is important to note that the definition 
in previous cosmological topology studies slightly differs from this. 
The genus, g, in these studies has been defined as the number of 
holes minus the number of connected regions: g = G — c. 



Fig. 2. — Illustration of alpha shapes. For two different values of a (top: large a; bottom: small a), we show the relation between the 
2-dimensional point distribution, the value of o, and the resulting alpha shape. Around each point in the sample, we draw a circle of radius 
a. The outline of the corresponding Voronoi tessellation is indicated by the edges (left). All Delaunay simplices dual to the decomposition 
of the union of disks by the Voronoi polygons are shown in black (vertices and edges) and red (triangles) (center). The resulting alpha 
shape is shown on the right. 



where the Euler characteristic x(dMi) is the integrated 
Gaussian curvature of the surface 

1 [ dr 

Here R\ {x) and i?2 (x) are the principal radii of curvature 
at point X of the surface. The integral of the Gaussian 
curvature is invariant under continuous deformation of 
the surface: perhaps one of the most surprising results 
in differential geometry. 

According to the Eulcr-Poincarc Formula, the Euler 
characteristic of the manifold M is equal to the alternat- 
ing sum of its Betti numbers, 

X(M) =/3o(M) -/3i(M)+/32(M). (4) 

For practical circumstances the third Betti number van- 
ishes, /33(M) = 0. The Euler characteristic of the bound- 
ary, 9M, is directly proportional to that of the Euler 
chara cteristic of the 3-manifold, M, (jSeifert fc Threlfall 
[l93l p.223) 

x(aM) = 2x(M). (5) 

The combination of these equations for the Euler char- 
acteristic establishes a fundamental relationship between 
differential geometry and algebraic topology. 

In the analysis described in this paper, we restrict 
ourselves to the three Betti numbers, /3o, /3i, and /32, 
of the 3-manifolds with boundary defined by the cos- 
mic mass distribution. In an upcoming paper, we study 
genus and Betti n umber characteristics of Gaussian fields 
(|Park et al.llMll) . 



3. ALPHA SHAPES 

One of the key concepts in the field of Computational 
Topology is alpha s hapes, introduced by Edelsbrun- 
ner and collaborators fEdclsbrunncr ct al."1983: 'Muccki 
[T993; Edclsbrunncr & Mucckc 1994; Edclsbrunncr 200f|. 
They generalize the convex hull of a point set and are 
concrete geometric objects that are uniquely defined for 
a particular point set S and a scale a. For their defini- 
tion, we look at the union of balls of radius a centered 
on the points in the set, and its decomposition by the 
corresponding Voronoi tessellation. The alpha complex 
consists of all Delaunay simplices that record the subsets 
of Voronoi cells that have a non-empty common intersec- 
tion within this union of balls (fig. 12. ip . 

The alpha complex thus consists of all simplices in the 
Delaunay triangulation that have an empty circumsphere 
with radius less than or equal to a. Here "empty" means 
that the open ball bounded by the sphere does not in- 
clude any points of S. For the extreme value a = 0, the 
alpha complex merely consists of the vertices of the point 
set. There is also a smallest value, amax, such that for 

^ Q!inax, the alpha complex is the Delaunay triangu- 
lation. The alpha shape is the union of simplices in the 
alpha complex and for a > amax the alpha shape is the 
convex hull of the point set. 

The alpha shape is a polytope in a fairly general sense: 
it can be concave and even disconnected. Its compo- 
nents can be three-dimensional clumps of tetrahedra, 
two-dimensional patches of triangles, one-dimensional 
strings of edges, and collections of isolated points, as well 
as combinations of these four types. 

Alpha shapes reflect the topological structure of a 
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Fig. 3. — Top: a shape of a GIF ACDM simulation, slice through 
simulation box. Bottom: a shape at 3 different values of a. 

point distribution on a scale parameterized by the real 
number a. Figure [3] shows the a-shape applied to a slice 
of a cosmological model. The three inserts show a-shapes 
at three different a values. 



3.1. Filtrations 

The concept of "filtration" is an important source of 
information about the topology of a point distribution. A 
filtration provides a view of the topology as a function of 
scale. Formally, given a space M, a filtration is a nested 
sequence of subspaces: 

= Mo C Ml C . . . C M,„ = M. (6) 

The nature of the filtrations depends, amongst other 
things, on the representation of the mass distribution. 
For a discrete point distribution, the set of alpha shapes, 
ordered by scale, constitute a filtration of the Delaunay 
tessellation. 

The set of real numbers a leads to a family of shapes 
capturing the intuitive notion of the overall versus fine 
shape of a point set. Starting from the convex hull grad- 
ually decreasing a, the shape of the point set gradually 
shrinks and starts to develop enclosed voids. These voids 
may join to form tunnels and larger voids. For negative 
a, the alpha shape is empty. 

Although the alpha shape is defined for all real num- 
bers a, there are only a finite number of different alpha 
shapes for any finite point set. In other words, the al- 
pha shape process proceeds discretely with increasing a, 
marked by the addition of new Delaunay simplices once 
a exceeds the corresponding threshold. 



Computing Betti Numbers 

algorithm in 2-D and 3-D for cal- 

of an a-complex is due to 

lDelfinad"o^ Edelsbrunner. (,1993< ). The algorithm is ef- 



3.2. 

The standard 
culating B etti numbers 



ficient in the sense that only one traversal of the sim- 
plices making up the a complex is required. First, the 
Delaunay tessellation of the point set is constructed. Its 
simplices - vertices, edges, triangles - can be ordered by 
the radius of the smallest sphere that touches the points 
defining the given simplex and contains no other data 
points. From this ordering a nested sequence of subcom- 
plexes can be constructed, each member of the sequence 
being constructed from the preceding one by the addition 
of a new simplex having a larger scale. If for example 
the added simplex is an edge, it will either connect two 
previously disjoint components or close a loop (thereby 
creating a tunnel). 

Cycling through all its simplices, in three dimensions 
the calculation is based on the following straightforward 
considerations. When a vertex is added to the alpha com- 
plex, a new component is created and (3q is increased by 
1. Similarly, if an edge is added, it connects two vertices, 
which either belong to the same or to different compo- 
nents of the current complex. In the former case, the 
edge creates a new tunnel, so /3i is increased by 1. In 
the latter case, two components get connected into one, 
so Po is decreased by 1. If a triangle is added, it either 
completes a void or it closes a tunnel. In the former 
case, f32 is increased by 1, and in the latter case, I3i is 
decreased by 1. Finally, when a tetrahedron is added, a 
void is filled, so P2 is lowered by 1. Following this pro- 
cedure, the algorithm has to include a technique for de- 
termining whether a p-simplex belongs to a p-eycle. For 
vertices and tetrahedra, this is rather trivial. For edges 
and triangles, a more elaborate procedure is necessary 
()Delfinado fc Edelsbrun^fTOOl . 

For software that implements these algorithmic ideas, 
we resort to the Computational Geometry Algorithms 
Libra ry, CGAL (|Caroh fc Teillaud,201k ^Da fc Yvine3 
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Fig. 4. — Evolution of the dark energy equation of state param- 
eter w(z) for ACDM (solid), SUGRA (dotted) and RP (dashed). 



4. TOPOLOGICAL ANALYSIS OF SIMULATIONS 
4.1. Our simulations 

CGAL is a C++ library of algorithms and data structures for 
computational geometry, see www . cgal . org, 
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Wc investigate the evolved structure in the standard 
ACDM model, the Ratra-Peebles (RP) model, and a 
SUGRA model. All simulations start with the primor- 
dial perturbation amplitudes and phases. The imprint 
of dark energy on the mass distribution is expressed via 
the different expansion history of the universe as a result 
of the differently evolving equation of state parameter 
w{z) (fig. |4]). This is determined by the evolution of the 
quintessence scalar field 0, forced by a parameterized 
potential Vjcj)). The DE model parameters (table 1 in 
IBos et all mm ) are chosen such that the models differ 
significantly from ACDM, while still in accordance with 
observations. For example, the slope of the quintessence 
potentials is set to ^rp = 0.347 and Csugra = 2.259. 

The models we consider are described by WMAP7 
(jLarson et al.ll201l[) best-fit value parameters (h = 0.704, 
rj„,. = 0.272, = 0.728, as = 0.809 and n, = 0.963). 
For each model we ran 512"^ and 256'^ particle simula- 
tions, in a periodic box of 200/i~^Mpc width, starting at 
z « 63. The baryonic content of the universe is incorpo- 
rated as dark matter particles, using flm = ^dm + = 
0.272. 

We use the SUBFIND algorithm (jSpringel et al.|[200l 
to identify the gravitationally bound haloes in the re- 
sulting particle distributions. At z = the simula- 
tion boxes contain 170740 (ACDM), 172949 (RP) and 
175332 (SUGRA) halos with mass higher than M - 
1.4 X 10"/i"1Mq (> 32 particles). They trace the gen- 
eral structures present in the field, even though a lot of 
the detailed structure is lost or diluted (fig. [T]). 

Slices from two simulations, at z = 0, are shown in fig- 
ure [TJ The left hand panel displays the standard ACDM 
model and the right hand panel shows the SUGRA 
model. The slices as shown in the upper panels are 
perceptually indistinguishable. Differences can only be 
perceived at higher resolution and lower contrast: this is 
shown in the lower panels of figure [TJ for the dark mat- 
ter density field as well as the haloes. Different filaments 
appear, the subclustering is different in the big cluster 
and there are a number of different voids. The voids dif- 
fer in size and shape, and thei r degree of merging (cf. 
ISheth fc van de Wevgaeri][200l . 

4.2. The analysis - Betti signatures 

We analyse the clustering of the haloes in these three 
simulations using Betti Numbers. In order to compute 
the Betti numbers for a point distribution in 3 dimen- 
sions, we generate the Delaunay tessellation from the 
point sample and from this the nested sequence of al- 
pha shapes. Subsequently, we determine the Betti num- 
bers of the set of alpha shapes as a function of a. The 
resulting Betti number versus a curves are the Betti sig- 
nature of the point distribution. In Figure 4 we show 
the Betti signatures for the scale dependence of /3o(a), 
/3i(a), and /32(q;) in the three simulations. The error bars 
in all panels are derived on the basis of estimates from 8 
sub-cubes. There are significant differences between the 
curves for /3i(a), and /32(q^) in the three models, show- 
ing that the Betti signature is a sensitive discriminator 
of structure, even in a single redshift slice. This remains 
true when taking into account differences in /3o and (Tg. 
In subsequent papers we follow the evolution of the Betti 
signature with time and assess their dependence on halo 





Fig. 5. — Betti signatures: Betti numbers, per unit volume, as 
a function of scale a at 2 = 0. Red: ACDM, Green: RP, Blue: 
SUGRA. Top: /3o; Middle: /3i; Bottom: ^2- Inserts: detailed 
differences near distribution peak. Lower strip: /3i//3i(ACDM). 
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mass and show that it provides a sensitive measure of 
the nature of dark energy. 

5. DISCUSSION AND CONCLUSIONS 

We have presented a powerful topology-based method 
for classifying the structure that develops in cosmologieal 
N-body simulations. The Betti signatures are an effective 
way of discriminating between models that superficially 
look the same. The method is an important generalisa- 
tion of the genus and Minkowski functional approaches, 
in that Betti analysis of alpha shapes is ab le to discrim- 
inate structural content at a deep er level (jRobinsI [20061 : 
Ivan de Wevgaert et al.llMollMll ). 

The Betti signatures of the cosmologieal simulations 
are shown to be capable of discriminating between mod- 
els that are difficult to distinguish by other means. This 
may involve dark energy signatures or traces of primor- 
dial non-Gaussianities in the cosmic mass distribution. 
This is because dark energy affects the expansion of the 
universe, which affects the formation of nonlinear ob- 
jects, while non-Gaussian initial conditions will change 
the formation history of galaxies and thus Betti num- 
bers. 

It is also important to put our approach in terms of 
Betti analysis in perspective relative to other measures 
of structure, in particular the Minkowski fun ctionals 
(|Mecke et aLlHOOl ISchmalzing fc Buchertlll997D . While 
the Euler characteristic and Betti numbers give informa- 
tion about the connectivity of a manifold, the other three 
Minkowski functionals are sensitive to local manifold de- 
formations. The Minkowski functionals give information 
about the metric properties of a manifold on the basis 
of which broad-brus h structural chara cteristics such as 
shape can be defined (jSahni et al.lll998D . The Betti num- 
bers focus on its topological properties, the scale depen- 
dence likewise providing some metric information. 

In a series of following papers we shall be 
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adding the powerful not ion of topological persistence 
( Edelsbrunner et aLl 12002} ) to our analysis of filtrations 
( Pranav et al.ll20l'lF . Betti analysis is an integral part of 
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1201 Ih . In standard practice, the multiscale nature of 
the mass distribution tends to be investigated by means 
of user-imposed filtering. Persistence entails the con- 
ceptual framework and language for separating scales of 
a spatial structure, and rationalizes the multiscale ap- 
proach by considering the range of filters at once. At 
the same time, it deepens the approach by combining 
it with topological measurements. Within the context 
of hierarchical cosmic structure formation, persistence 
therefore provides a natural formalism for a multiscale 
topology study of the Cosmic Web. Alpha shapes pro- 
vide the perfect context for understanding the concept: 
the position of a feature within the structural hierarchy is 
determined by the a interval over which it persists. Per- 
sistence even provides a natural path towards removing 
topological noise, which would be identified as feature s 
with small persistence (see Ivan de Wevgaert et al.ll201ll) . 
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